CETA  80-5 


Interpretation  of  Wave  Energy  Spectra 


Edword  F.  Thompson 


COASTAL  ENGINEERING  TECHNICAL  AID  NO.  80-5 

JULY  1980 


IS  | j  Sj 


Approved  for  public  release; 
distribution  unlimited. 


DT1C 

.ELECTE 


SEP  2  9  1980 


U.S.  ARMY,  CORPS  OF  ENGINEERS 

COASTAL  ENGINEERING 
RESEARCH  CENTER 

Kingman  Building 
Fort  Belvoir,  Va.  22060 


80  9  26  02 


Reprint  or  republication  of  any  of  this  material  shall  give  appropriate 
credit  to  the  U.S.  Army  Coastal  Engineering  Research  Center. 

Limited  free  distribution  within  the  United  States  of  single  copies  of 
this  publication  has  been  made  by  this  Center.  Additional  copies  are 
available  from: 

National  Technical  Information  Service 
A  TTN:  Operations  Division 
5285  Port  Royal  Road 
Springfield,  Virginia  22161 

The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated  by  other 
authorized  documents. 


on. II.. _ _ 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  PM.  Entered) 

REPORT  DOCUMENTATION  PAGE  befoIe^^etogtorm 

I.  REPORT  NUMBER  I*.  GOVT  ACCESSION  NO.  »•  RECIPIENT’S  CATALOG  NUMBER 


CETA  80-5 

4.  TITLE  (tnd  Submit) 


ifjP' 


ERPRETATION  OF  WAVE  ENERGY  SPECTRA 


dk  Coastal 
.  *L)  Technical 


r  WErPRT  »  PFHirtn  rny^mwi 

al  Engineering  , 
ical  Aid  PI* 


17.  AUTHORfr) 


I.  CONTRACT  OR  GRANT  HUMBERT*.) 


Edward  F .  /Thompson 


7 


S.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS  10.  PROGRAM  ELEMENT.  PROJECT.  TASK 

Department  of  the  Army  *REA  *  W0Rlt  0H,T  bombers 

Coastal  Engineering  Research  Center  (CERRE-CO)  A31463 

Kingman  Building,  Fort  Belvoir,  Virginia  22060 

II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS  - REPORT-DRY* - - 

Department  of  the  Army  f  fl )  JulfH§||8j/  / 

Coastal  Engineering  Research  Center  *•>  number  o»-p abes-^3 

Kingman  Building,  Fort  Belvoir,  Virginia  22060 _ 21 _ 

<4.  MONITORING  AGENCY  NAME  4  AOORESy/l  dMtrent  ham  Controlling  OfUct)  IS.  SECURITY  CLASS,  (at  »'•  upon) 


Jul*fl§8G 


JljZBRt-tETA'  It'* 


UNCLASSIFIED 


19«.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


1 16.  DISTRIBUTION  STATEMENT  ( oi  thim  Rmport) 


Approved  for  public  release;  distribution  unlimited. 


I  17.  DISTRIBUTION  STATEMENT  (Ol  the  tbttrtcl  entered  In  Block  20,  It  tfUbMlI  Bom  Report) 


IIS.  SUPPLEMENTARY  NOTES 


I  ts.  KEY  WORDS  (Continue  on  rtrtrtt  tide  It  mniw jr  md  Identify  br  block  mmber) 


Spectral  analysis 
Wave  energy  spectra 


Wave  gages 
Wave  height 


Wave  hindcast 
Wave  period 


to.  ^DETRACT  fCKM»  cm  root no  N*  R  Hoeooeen  mi  Idemtltr  bybleek  lumber) 

^Guidelines  for  interpreting  nondirectional  wave  energy  spectra  are  presented. 
A  simple  method  is  given  for  using  the  spectrum  to  estimate  a  significant  height 
and  period  for  each  major  wave  train  in  most  sea  states.  The  method  allows  a 
more  detailed  and  accurate  description  of  ocean  surface  waves  than  that  given  by 
a  single  significant  height  and  period,  yet  it  eliminates  much  of  the  formidable 
detail  of  a  full  spectrum.  An  example  problem  illustrating  application  of  the 
method  is  presented.  Spectral  analysis  and  display  techniques,  and  the  natural 
variation  of  spectra  in  space  and  time,  are  discussed. _ _ _ _ 


no  « 

w  i  ji 


EDITION  OP  I  NOV  ••  IS  I 


5SIFIED 


SECURITY  CLMStrtCAT 


.$6 


PREFACE 


This  report  presents  guidelines  for  interpreting  nondirectional  wave  energy 
spectra.  The  guidelines  apply  to  spectra  derived  from  both  wave  gage  measure¬ 
ments  and  from  numerical  wave  hindcasting  models.  A  method  is  provided  for 
using  the  spectrum  to  estimate  a  significant  height  and  period  for  each  major 
wave  train  in  a  sea  state,  except  major  wave  trains  with  nearly  the  same  period 
and  different  directions  cannot  be  distinguished.  The  method  has  undergone 
limited  testing  and  has  been  applied  to  7  station-years  of  gage  data,  but  fur¬ 
ther  testing  in  we 11 -documented  field  situations  is  needed.  The  guidelines  and 
method  are  consistent  with  but  are  more  practical  and  explicit  than  the  material 
in  the  Shore  Protection  Manual  (SPM) .  The  work  was  done  under  the  wave  measure¬ 
ment  program  of  the  U.S.  Army  Coastal  Engineering  Research  Center  (CERC) . 

This  report  was  prepared  by  Edward  F.  Thompson,  Hydraulic  Engineer,  under 
the  general  supervision  of  Dr.  C.L.  Vincent,  Chief,  Coastal  Oceanography  Branch. 
Helpful  reviews  by  Dr.  C.L.  Vincent,  Dr.  D.L.  Harris,  P.  Knutson,  and  P.  Vitale 
are  acknowledged.  Dr.  D.  Esteva  provided  the  data  for  Figure  6. 

Comments  on  this  publication  are  invited. 


Approved  for  publication  in  accordance  with  Public  Law  1966,  79th  Congress, 
approved  31  July  1945,  as  supplemented  by  Public  Law  1972,  88th  Congress, 
approved  7  November  1963. 


Colonel,  Corps  of  Engineers 
Commander  and  Director 
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CONVERSION  FACTORS ,  U.S.  CUSTOMARY  TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 


U.S.  customary  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  (SI)  units  as  follows: 


Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x  10-3 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton ,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.01745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

*To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  formula:  C  ■  (5/9)  (F  -32). 


To  obtain  Kelvin  (K)  readings,  use  formula:  K  ■  (5/9)  (F  -32)  ♦  273.15. 


SYMBOLS  AND  DEFINITIONS 


d 

E 


Ej »  E(fj) 


L-mor 


rpi 


H 


'St 


wave  component  amplitude 
water  depth 
sum  of  all  Ej 

spectral  energy  density  values 

energy  density  for  the  highest  spectral  peak 

frequency  of  spectral  component  in  hertz 

frequency  corresponding  to  the  highest  spectral  peak 

frequency  corresponding  to  the  ith  highest  spectral  peak 

acceleration  due  to  gravity 

significant  wave  height  corresponding  to  the  full  spectrum 

significant  wave  height  corresponding  to  the  ith  highest  spectral 
peak 


Kl,  K2 


K1 K2 i 


N 


*ty 

Spi 


Af,  (Af)j* 


n 


indices  representing  the  lower  and  upper  bounds  of  the  highest 
spectral  peak 

indices  representing  the  lower  and  upper  bounds  of  the  ith  highest 
spectral  peak 

total  number  of  spectral  frequency  components 
spectral  energy  values 

energy  contained  in  ith  highest  spectral  peak 
time 

period  corresponding  to  the  highest  spectral  peak 

period  corresponding  to  the  ith  highest  spectral  peak 

frequency  bandwidth  represented  by  each  spectral  energy  density 
(in  hertz) 

sea-surface  elevation  referenced  to  local  mean  water  level 
phase  of  spectral  component 

frequency  of  spectral  component  (in  radians  per  second) 
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INTERPRETATION  OF  HAVE  ENERGY  SPECTRA 


by 

Eduard  F.  Thompson 
I  INTRODUCTION 

The  ocean  usually  has  more  than  one  independent  train  of  waves  propagating 
along  its  surface  in  U.S.  coastal  areas.  The  common  practice  of  using  a  single 
significant  height  and  period  for  a  sea  state  can  be  misleading  because  no  indi¬ 
cation  is  given  to  the  existence  or  characteristics  of  other  trains.  On  the 
other  hand,  an  estimate  of  the  wave  energy  spectrum  provides  more  information 
than  is  generally  used  in  coastal  engineering.  The  spectrum  can  be  reduced  to 
estimates  of  significant  wave  height  and  period  for  all  major  wave  trains  pres¬ 
ent.  A  knowledge  of  these  characteristics  for  major  wave  trains  is  often  im¬ 
portant  to  coastal  engineers. 

Spectra  are  becoming  widely  available  through  various  field  wave  measurement 
programs,  laboratory  tests  with  programable  wave  generators,  and  numerical  wave 
hindcasting  projects.  Because  of  the  availability  and  applications  of  spectra, 
practicing  coastal  engineers  should  become  familiar  with  spectra  and  their 
interpretation. 


II.  PRIMARY  GOAL  OF  SPECTRAL  ANALYSIS 

A  fundamental  parameter  for  characterizing  a  wave  field  is  some  measure  of 
the  periodicity  of  the  waves.  For  many  years  a  significant  period,  which  could 
be  subjectively  estimated  in  various  ways,  was  used.  However,  the  ocean  surface 
often  has  waves  characterized  by  several  distinct  periods  occurring  simultane¬ 
ously.  A  record  of  the  variation  of  sea-surface  elevation  with  time,  commonly 
called  a  time  series,  frequently  appears  confusing  and  is  difficult  to  interpret. 

Developments  in  computer  technology  and  in  mathematical  analysis  of  time 
series  have  provided  a  practical  approach  to  an  objective,  more  comprehensive 
analysis  of  periodicity  in  wave  records.  The  approach  is  to  express  the  time 
series  as  a  sum  of  periodic  functions  with  different  frequency,  amplitude,  and 
phase.  The  simplest  functions  to  apply  are  the  trigonometric  sine  and  cosine 
functions.  Thus,  the  time  series  of  sea-surface  deviations  from  the  mean  sur¬ 
face,  n(t),  is  expressed  by  equation  (3-11)  in  the  Shore  Protection  Manual 
(SPM)  (U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering  Research  Center,  1977) 
as 

nCt)  -  ^  aj-  cos  (mjt-ty)  (1) 

where 

a j  =  amplitude 

u j  =  frequency  in  radians  per  second 
*  phase 


Frequency  is  often  expressed  in  terms  of  hertz  units  where  one  hertz  is  equal 
to  one  cycle  per  second.  One  hertz  is  also  equivalent  to  2ir  radians  per 
second.  If  the  symbol  denotes  frequency  in  hertz,  then  2nf^  *  uj. 

The  amplitudes,  aj,  computed  for  a  time  series,  give  an  indication  of  im¬ 
portance  of  each  frequency,  fj*.  The  sum  of  the  squared  amplitudes  is  related 
to  the  variance  of  sea-surface  elevations  in  the  original  time  series  (eq.  3-12 
in  the  SPM)  and  hence  to  the  potential  energy  contained  in  the  wavy  sea  surface. 
Because  of  this  relationship,  the  distribution  of  squared  amplitudes  as  a  func¬ 
tion  of  frequency  can  be  used  to  estimate  the  distribution  of  wave  energy  as  a 
function  of  frequency.  This  distribution  is  called  the  energy  spectrum  and  is 
often  expressed  as 

a2. 

(Ej)mj  =  4=  sj  (2) 

where  Ej  =  E  (fj)  is  the  energy  density  in  jth  component  of  energy  spectrum, 
(Af)j'  the  frequency  bandwidth  in  hertz  (difference  between  successive  fj) , 
and  Sj  =  S (fj)  energy  in  jth  component  of  energy  spectrum.  Equation  (2)  is 
similar  to  equation  (3-15)  in  the  SPM. 


An  energy  spectrum  computed  from  an  ocean  wave  record  is  plotted  in  Fig¬ 
ure  1.  Frequencies  associated  with  large  values  of  energy  density  (or  large 
values  of  a?/[2(Af)^]l  see  eq.  2)  represent  dominant  periodicities  in  the  orig¬ 
inal  time  series.  Frequencies  associated  with  small  values  of  a^/[2(Af)  •]  are 

d  d 

usually  unimportant.  It  is  common  for  ocean  wave  spectra  to  show  two  or  more 
dominant  periodicities  (Fig.  1). 


jj  15,000 


<=>  10,000 


fp*  fp  Frequency  (Hz) 

Figure  1.  Spectrum  for  Wrightsville  Beach,  North  Carolina,  0700  e.s.t., 

12  February  1972;  H*  ■  4.2  feet  (128  centimeters),  Af  -  0.01074 
hertz,  and  depth  ■  17.7  feet  (5.4  meters). 


The  primary  goal  of  spectral  analysis  is  to  objectively  identify  all  im¬ 
portant  frequencies  in  a  wave  record.  Since  wave  period  in  seconds  is  equal 
to  the  reciprocal  of  frequency  in  hertz,  important  wave  periods  are  also 
identified. 

III.  PRACTICAL  LIMITATIONS  OF  SPECTRAL  ANALYSIS 


1.  Calculation  Procedures. 

The  appearance  of  a  spectrum  can  be  noticeably  influenced  by  the  methods 
used  for  calculation  and  display,  neither  of  which  is  standardized  in  coastal 
and  ocean  engineering  activities. 

Spectra  are  computed  from  both  digital  wave  records  and  analog  records  for 
which  an  assortment  of  analog  spectral  analysis  devices  exists.  A  spectrum  is 
computed  from  a  digital  wave  record  by  either  (a)  computing  the  fourier  trans¬ 
form  of  the  autocovariance  function  of  the  record,  or  (b)  computing  the  fourier 
transform  of  the  record  directly  from  the  record  using  the  fast  fourier  trans¬ 
form  (FFT)  approach  (see  Harris,  1974,  for  further  detail).  Both  of  these 
algorithms  are  used  in  conjunction  with  wave  records  subjected  to  various  fil¬ 
ters  and  smoothing  functions  before  analysis.  Further,  some  form  of  smoothing, 
averaging,  or  summing  is  often  applied  to  the  computed  spectral  components. 

Different  methods  for  calculating  a  spectrum  will  produce  slightly  different 
estimates  of  the  spectrum  when  applied  to  a  particular  wave  record.  Major  dif¬ 
ferences  in  the  height  of  the  spectral  peak  were  shown  by  Wilson,  Chakrabarti, 
and  Snider  (1974)  when  different  approximations  to  the  autocovariance  function 
were  used  and  different  smoothing  functions  were  applied  to  the  spectrum  of  a 
field  wave  record  (Fig.  2).  Major  differences  in  the  height  of  the  spectral 
peak  and  energy  levels  between  peaks  were  noted  by  Harris  (1974)  when  a  time 
series  composed  of  three  superimposed  sinusoidal  waves  was  analyzed  by  several 
accepted  methods  (Fig.  3). 

Hindcast  wave  energy  spectra  are  confuted  by  estimating  atmospheric  input 
of  energy  to  the  sea  surface  and  redistribution  of  energy  within  a  spectrum. 

The  estimates  are  based  on  a  series  of  equations  derived  from  the  physics  of 
air-sea  interaction  and  waves.  The  quality  and  characteristics  of  hindcast 
spectra  are  a  function  of  the  model  used  to  perform  the  calculations  (Resio 
and  Vincent,  1979)  as  well  as  the  accuracy  of  the  input  wind  field. 

Spectra  obtained  from  either  measurements  or  hindcasts  are  also  limited  by 
the  resolution  of  the  computation  technique.  The  energy  density  and  frequency 
at  a  spectral  peak  can  be  noticeably  distorted  if  the  frequency  bandwidths, 
(Af)j,  are  not  small  enough  to  permit  clear  definition  of  major  peaks. 

2.  Display  Formats. 

The  appearance  of  a  spectrum  can  be  strongly  affected  by  the  display  format 
used.  Harris  (1972)  showed  five  often  used  formats  for  plotting  spectra  (Fig. 

4) .  Each  format  alters  the  appearance  of  the  spectrum.  Format  E  shifts  the 
relative  magnitudes  of  spectral  peaks  enough  that  the  second  highest  peak  in 
formats  A,  B,  C,  and  D  becomes  the  highest  peak  in  format  E.  A  and  C  are  the 
two  most  frequently  used  formats. 


Frequency  (  Hi  ) 


Figure  2 


Figure  3 


Energy  spectra  for  wave  record  at  0400  e.s.t.,  17  March 
1968,  from  a  weather  ship  in  the  North  Atlantic,  computed 
with  different  approximations  to  the  autocovariance  func¬ 
tion  and  different  spectral  smoothing  functions  (after 
Wilson,  Chakrabarti,  and  Snider,  1974). 


Energy  spectra  for  record  composed  of  three  superimposed 
sinusoidal  waves.  Simulated  frequencies  indicated  by 
vertical  lines  from  top  of  graph  (after  Harris  1974). 


Figure  4.  Five  formats  frequently  used  in  displaying  wave  energy  spectra.  The 
actual  spectrum  is  identical  in  all  five  graphs.  The  frequency  band 
width,  (Af)j,  is  constant  for  all  j  so  that  energy  values  and  energy 
density  values  differ  by  a  constant  factor  (see  eq.  2).  The  plots 
would  look  the  same  if  energy  were  replaced  by  energy  density,  but 
the  vertical  scale  would  change  (from  Harris,  1972). 


3.  Natural  Variability. 

Wave  energy  spectra  are  naturally  variable  simply  because  they  are  based  on 
a  finite  length  record  of  a  wave  field  which  varies  in  time  and  space.  Spectra 
computed  for  successive  records  of  a  relatively  stationary  wave  field  are  never 
identical  and  often  differ  noticeably.  The  magnitude  of  spectral  variation  in 
time  is  illustrated  by  spectra  derived  at  2-hour  intervals  from  two  pressure 
gages  along  the  southern  California  coast  (Fig.  5).  The  significant  wave  height 
is  nearly  constant  in  the  figure. 


Spatial  variation  of  the  spectrum  over  short  alongshore  distances  in  shallow 
water  is  also  shown  in  Figure  5.  Each  spectrum  in  the  top  row  of  the  figure  can 
be  compared  to  the  spectrum  immediately  below  it  to  see  variations  between  spec¬ 
tra  from  two  gages  80  feet  (24  meters)  apart.  In  this  figure,  spatial  variations 
are  smaller  than  temporal  variations.  Spatial  variations  would  be  expected  to 
be  greater  if  the  gages  were  farther  apart  or  the  water  depth  varied  between 
measurement  points.  Variations  between  spectra  from  gages  situated  along  a  line 
perpendicular  to  shore  are  shown  in  Figure  6.  The  spatial  variations  are  more 
prominent  in  this  figure  than  in  Figure  5. 


II 


Entrgy  Oensity  Entrgy  Otnsity 


Frtqutncy  (  Hi)  Frtqwncy  (  Hi ) 

Figure  6.  Wave  energy  spectra  from  pier-mounted  continuous  wire  gages 
at  the  CERC  Field  Research  Facility  (FRF)  near  Duck,  North 
Carolina,  showing  variation  along  a  line  perpendicular  to 
shore.  Solid  lines  represent  a  gage  at  the  seaward  pier  end 
(depth  29  feet  or  8.8  meters);  dashlines  represent  a  gage  480 
feet  (146  meters)  from  the  seaward  pier  end  (depth  22  feet 
or  6.6  meters);  dot-dash  lines  represent  a  gage  840  feet  (256 
meters)  from  the  pier  end  (depth  17  feet  or  5.1  meters). 
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IV.  INTERPRETATION  OF  SPECTRA 


Because  of  natural  variability  in  the  spectrum  and  artificial  variability 
induced  by  analysis  and  display  techniques ,  the  spectrum  should  never  be  in¬ 
terpreted  as  an  exact  representation  of  energy  density  versus  frequency  for  a 
wave  field.  However,  certain  major  features  of  the  spectrum  are  consistent 
and  meaningful. 

1 .  Highest  Spectral  Peak. 

a.  Frequency  and  Period.  Frequency  corresponding  to  the  highest  spectral 
peak,  fp,  is  usually  a  reliable  measure  of  the  dominant  wave  frequency;  fp 
is  shown  in  Figure  1.  Period  corresponding  to  the  highest  spectral  peak.  Ip, 
is  equal  to  the  reciprocal  of  fp  and  is  usually  a  good  estimate  of  the  domi¬ 
nant  wave  period. 

b.  Energy  and  Significant  Wave  Height.  Energy  contained  in  the  highest 
peak,  Spi,  is  defined  as  the  total  energy  in  the  vicinity  of  the  highest  peak. 

S2 

Spl  »  l  E^flj  (3) 

j=Ki+l 


where  Kj  and  K2  are  indices  representing  the  lower  and  upper  bounds  of  the 
main  peak.  The  upper  and  lower  bounds  sometimes  represent  a  broad  range  of 
frequencies  (see  Fig.  1).  Spl  is  relatively  consistent,  and  is  less  influenced 
than  the  magnitude  of  the  highest  peak  by  data  collection  procedures,  by  analy¬ 
sis  and  summarization  procedures,  and  by  temporal  and  spatial  variation. 


Some  spectral  analysis  procedures  are  designed  so  that  (Af  ) j 
constant  for  all  j,  which  leads  to 

*2 

Spl  =  (Af)  l  Ej 
j=Ki+l 


Af  is  a 


(3a) 


Significant  wave  height  corresponding  to  highest  spectral  peak  Hai,  is  an 
estimate  of  the  significant  height  for  the  wave  train  represented  by  the  highest 
spectral  peak.  It  is  computed  by  the  relationship 

Hsl  *  4V"Spl  (4) 

Energy  density  at  the  highest  spectral  peak,  E max*  can  be  an  indicator  of 
how  well  focused  the  wave  energy  is  in  frequency.  Although  this  parameter  is 
variable,  major  differences  in  E max  (°n  the  order  of  SO  percent)  between 
spectra  analyzed  by  the  same  method  can  be  meaningful.  Emax  is  shown  in 
Figure  1. 


2.  Major  Secondary  Spectral  Peaks. 

a.  Identification  of  Major  Secondary  Peaks.  Major  secondary  spectral  peaks 
are  often  indicative  of1  independent  secondary  wave  trains  characterized  by  dif¬ 
ferent  heights,  periods,  and  directions  than  the  train  represented  by  the  main 
peak  (examples  are  given  in  McClenan  and  Harris,  1975).  Identification  of  major 
secondary  spectral  peaks  involves  some  subjective  judgment,  but  an  objective 
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test  for  major  secondary  peaks  has  been  developed  and  used  at  CERC.  The  test 
is  applied  to  the  difference  in  energy  density  between  a  spectral  peak  and  the 
preceding  and  following  spectral  valleys.  If  that  difference  exceeds  3  percent 
of  the  total  of  all  spectral  energy  density  values,  E,  then  the  peak  is  con¬ 
sidered  major.  Details  of  the  procedure  with  a  computerized  version  are  given 
in  Thompson  (1980).  The  procedure  was  applied  to  the  spectrum  in  Figure  1,  com 
puted  from  an  ocean  wave  record,  and  two  major  peaks  were  identified.  The  pro¬ 
cedure  has  been  applied  to  7  station-years  of  shallow-water  spectra  by  Thompson 
(1980)  to  show  that  two-thirds  of  the  ocean  and  gulf  coast  spectra  have  more 
than  one  major  peak. 

If  two  independent  trains  have  nearly  the  same  frequency  but  different  di¬ 
rections,  they  cannot  be  identified  by  the  method  presented  in  this  report. 
However,  an  analogous  method  could  be  developed  for  use  with  directional  wave 
spectra  to  identify  all  wave  trains  which  have  distinct  frequencies  or  direc¬ 
tions.  Directional  spectra  which  give  estimates  of  the  distribution  of  wave 
energy  density  as  a  function  of  both  frequency  and  direction  are  becoming  in¬ 
creasingly  available,  mainly  from  spectral  hindcasting  models  but  also  from 
improved  measurement  devices.  A  directional  spectrum  obtained  by  the  Data  Buoy 
Office  of  the  National  Oceanic  and  Atmospheric  Administration  (NOAA)  from  an 
experimental  large  discus  buoy  (Burdette,  Steele,  and  Trampus,  1978)  is  shown 
in  Figure  7.  The  figure  indicates  a  concentration  of  low-frequency  wave  energy 
coming  from  the  north-northeast  and  a  concentration  of  high-frequency  energy 
coming  from  the  quadrant  between  east  and  south. 


Figure  7.  Directional  spectrum  obtained  in  the  Atlantic  Ocean  68  miles  (110 
kilometers)  east  of  Jacksonville,  Florida,  at  1900  e.s.t.,  30  March 
1977.  Contour  values  are  in  units  of  0.001  meter  squared  per 
hertz  per  degTee.  The  radial  coordinate  is  frequency  in  hertz;  the 
azimuthal  coordinate  is  direction  from  which  energy  is  coming  (from 
Burdette,  Steele,  and  Trampus,  1978). 


Identification  of  major  spectral  peaks  can  also  be  a  useful  step  toward 
approximating  water  velocities  and  accelerations  from  a  complex  record  of  sea- 
surface  elevations,  although  procedures  for  doing  this  are  not  well  established. 
Velocities  and  accelerations  could  be  estimated  for  each  peak  separately,  using 
equations  in  the  SPM  (Sec.  2.234);  then  they  could  be  added  together  to  give 
estimates  of  resultant  velocity  and  acceleration.  If  direction  estimates  are 
available,  the  velocity  and  acceleration  components  should  be  added  vectorially. 
An  alternative  approach  is  to  use  the  nondirectional  spectrum  to  estimate  a 
velocity  spectrum  (Harris,  1972). 

Occasionally,  spectra  have  major  secondary  peaks  which  do  not  represent 
independent  wave  trains.  In  particular,  when  waves  are  very  steep  or  in  very 
shallow  water  the  spectrum  will  often  have  secondary  peaks  at  frequencies  which 
are  integral  multiples  of  the  dominant  frequency.  The  secondary  peaks  indicate 
that  additional  frequencies  are  needed  to  represent  the  nonsinusoidal  wave  pro¬ 
files  (the  peaks  do  not  necessarily  represent  independent  wave  trains).  When 
a  major  secondary  peak  appears  at  twice  the  dominant  frequency,  its  source  can 
usually  be  identified  by  referring  to  a  plot  of  the  time  series  from  which  the 
spectrum  was  computed.  The  time  series  is  expected  to  show  wave  profiles  which 
are  clearly  nonsinusoidal  if  the  major  secondary  peak  is  non independent.  Coarse 
guidelines  for  when  nonindependent  secondary  peaks  may  occur  are 

steepness:  >  0.008  (5) 

g Tp 

or 

relative  depth:  <  0.01  (6) 

gTp 


In  cases  where  both  steepness  and  relative  depth  approach  the  above  guidelines, 
non independent  peaks  may  also  be  evident. 

b.  Frequency  and  Period.  Frequency  and  period  corresponding  to  major 
secondary  spectral  peaks,  f^£  and  T usually  indicate  secondary  frequen¬ 
cies  and  periods  at  which  major  amounts  of  energy  are  present.  fp2,  corre¬ 
sponding  to  the  second  highest  spectral  peak,  is  shown  in  Figure  1.  Equations 
(5)  and  (6)  can  be  helpful  in  identifying  cases  where  major  secondary  peaks  do 
not  represent  independent  energy  concentrations. 

c.  Energy  and  Significant  Wave  Height.  Energy  and  significant  wave  height 
corresponding  to  major  secondary  spectral  peaks,  Spi,  and  Hgi,  indicate  the 
relative  importance  of  secondary  peaks.  Since  major  secondary  peaks  often  rep¬ 
resent  independent  wave  trains,  Spi  and  H si  can  be  estimates  of  the  energy 
and  significant  height  of  secondary  trains.  A  method  used  at  CERC  to  estimate 
S pi  is  to  partition  the  spectrum  at  its  lowest  point  between  every  pair  of  ad¬ 
jacent  major  peaks  (Fig.  1).  The  energy  between  partitions  is  then  totaled  by 


K2i 

h*  "  .  1  ,  (7) 

where  Kit  is  the  index  representing  the  lower  bound  of  the  peak,  and  K2t  the 
index  representing  the  upper  bound  of  the  peak.  If  (Af)j'  »  4f  »  constant,  then 
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(7a) 


Kz  i 

Spi  •  (Af)  l  Ej 


Significant  height  is  estimated  from  the  total  energy  assigned  to  each  peak  by 
an  equation  similar  to  equation  (4) . 

H ei  =  4^  (8) 

3.  Example  Problem. 

This  problem  illustrates  a  method  for  estimating  peak  frequencies,  periods, 
and  significant  heights  for  a  spectrum  with  two  major  peaks. 

GIVEN :  Wave  spectrum  from  Huntington  Beach,  California,  for  which  the  signifi¬ 
cant  height  is  Hs  =  5.7  feet  (175  centimeters)  (Fig.  8  and  Table).  Energy 
density  is  expressed  as  a  percent  of  the  sum  of  energy  densities  for  all  f j 
listed  in  Table.  Hs  =  5.7  feet  is  based  on  100  percent  of  the  energy  in  the 
spectrum. 

FIND:  Estimate  a  separate  significant  wave  height  and  peak  period  for  each 
wave  train  indicated  by  the  spectrum. 

SOLUTION:  To  identify  major  spectral  peaks,  the  difference  in  energy  density 
between  peak  2  and  the  valley  between  peaks  1  and  2  is  estimated  along  the 
vertical  axis  in  Figure  8  or  from  Table  (about  4.5  percent).  Since  this  is 
greater  than  3  percent,  peaks  1  and  2  are  accepted  as  major  peaks.  Several 
other  peaks  appear  at  frequencies  higher  than  0.15  hertz,  but  there  is  no 
other  combination  of  peak  and  valley  for  which  the  difference  in  energy  den¬ 
sity  exceeds  3  percent.  Therefore,  peaks  1  and  2  are  the  only  major  peaks. 
Frequencies  for  peaks  1  and  2  are  estimated  along  horizontal  axis  or  from 
Table : 


fp  «  f8  =  0.081  Hz 
ip 2  *=  f13  =  0.135  Hz 


Since  fp2  is  not  integral  multiple  of  fp,  assume  that  peaks  1  and  2  repre¬ 
sent  independent  wave  trains.  The  reciprocal  of  each  peak  frequency  gives 
peak  period. 


Tpz 


=»  _L  =  12.3 
fP 


s 


s 


Compute  total  energy  density,  E,  in  the  full  spectrum  by  combining  equations 
(3a)  and  (4)  and  rearranging  to  get 


E  = 


-  192  ft2/Hz 


Energy  Density  ( pet ) 


Peak  l 


1  1  Region  used  to  Compute  Spi 

Region  used  to  Compute  Sp? 


0.2  0.3 

Frequency  (Hz) 


0.4 


fp  fp2 

Figure  8.  Spectrum  for  Huntington  Beach,  California,  0640  e, 
24  August  1972;  H*  «  S.7  feet  (17S  centimeters), 
Af  ■  0.01074  hertz,  and  water  depth  =  26  feet  (2.t 


Table.  Spectrum  for  Huntington  Beach,  California, 
0640  e.s.t. 


i 

Td 

(Hz) 

(pet) 

j 

- rj 

<IU) 

EJ 

(pet) 

l 

0.005 

0 

25 

0.264 

0.4 

2 

0.016 

0 

26 

0.274 

0.3 

3 

0.027 

0 

27 

0.285 

0.3 

4 

0.038 

0.2 

28 

0.296 

0.4 

S 

0.049 

1.1 

29 

0.307 

0.4 

6 

0.059 

5.1 

30 

0.317 

0.2 

7 

0.070 

19.3 

31 

0.328 

0.1 

8 

0.081 

21.9 

32 

0.339 

0.3 

9 

0.092 

7.0 

33 

0.349 

0.4 

10 

0.102 

5.9 

34 

0.360 

0.1 

11 

0.113 

3.1 

35 

0.371 

0.2 

12 

0.124 

5.9 

36 

0.382 

0.1 

13 

0.135 

7.6 

37 

0.392 

0.3 

14 

0.14S 

4.9 

38 

0.403 

0.3  j 

IS 

0.156 

2.6 

39 

0.414 

0.2 

16 

0.167 

1.3 

40 

0.425 

0.1 

17 

0.178 

1.8 

41 

0.435 

0.4 

18 

0.188 

0.7 

42 

0.446 

0.3 

19 

0.199 

0.7 

43 

0.457 

0.2 

20 

0.210 

0.5 

44 

0.468 

0.1 

21 

0.221 

0.4 

45 

0.478 

0.2 

22 

0.231 

1.3 

46 

0.489 

0.1 

23 

0.242 

0.6 

47 

0.500 

0.1 

24 

0.253 

0.4 

48  to  100 

0. 5-1.0 

2.2 

*HS  «  5.7  feet  (175  ccntiwetcrs) ;  if  ■  0.01074  Hertz. 


0.5 

s.t. , 
meters) . 


Estimate  the  energy  in  peak  1  using  equation  (3a),  the  tabulated  energy  den¬ 
sity  (Table) ,  and  the  above  value  of  E  to  convert  the  sum  of  energy  density 
from  percent  to  foot  squared. 

11 

Spl  *  (Af)  l  Ej 
3*4 

=  (0.01074  Hz) 

=  1.31  ft2 

Similarly  estimate  the  energy  in  peak  2  using  equation  (7a) 

100 

Sp2  =  (Af)  l  Ej 
j=12 

=  (0.01074  Hz)  (3-^o  pc't)^192  ft2/Hz) 

=  0.75  ft2 

Estimate  significant  height  for  peak  1  using  equation  (4) . 

Hel  =  4  V®pl 

*  4 V 1-31  ft2  *  4.6  ft 

Note  that  Hgl  is  1.1  feet  lower  than  Hg  based  on  the  full  spectrum.  Simi¬ 
larly,  significant  height  for  peak  2  is  estimated  from  equation  (8) 

Hs2  *  3.5  ft 

Hs,  Hsi  and  Hs2  are  related  in  this  example  by 

H«  %/h2!  ♦  H22 

V.  INTERPRETATION  OF  SPECTRA  FOR  APPLICATIONS 
SENSITIVE  TO  SPECIFIC  FREQUENCIES 

The  interpretation  of  spectra  in  terms  of  wave  trains  is  appropriate  for 
most  coastal  engineering  work.  However,  certain  engineering  applications,  es¬ 
pecially  applications  in  which  resonance  can  occur,  are  highly  sensitive  to  one 
frequency  or  a  small  range  of  frequencies.  Estimates  of  how  much  energy  can  be 
expected  at  that  frequency  or  range  of  frequencies  are  required.  The  estimates 
are  obtained  directly  from  the  spectrum.  Since  the  estimates  are  sensitive  to 
data  collection  and  analysis  procedures,  it  is  especially  important  that  the 
procedures  be  optimum  when  such  applications  are  intended  or  anticipated. 

For  frequency-sensitive  applications,  it  is  generally  assumed  that  each 
frequency  with  nonzero  energy  represents  an  independent  wave  component,  regard¬ 
less  of  whether  it  is  a  spectral  peak.  Thus,  in  a  resonance  problem  sensitive 
to  frequency  f^,  spectral  energy  at  f#  is  treated  as  an  independent  wave. 

If  fp  corresponds  to  a  major  spectral  peak,  it  would  be  appropriate  (and 
conservative)  to  estimate  the  energy  from  equation  (7).  In  a  floating  break¬ 
water  problem  where  energy  at  frequencies  lower  than  some  cutoff  frequency, 


(t5T^)<192  “2/hz) 
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f<5,  will  be  transmitted  through  the  breakwater,  the  transmitted  energy  is  esti¬ 
mated  from  an  incident  spectrum  by  summing  all  energy  at  frequencies  lover  than 
fa,  regardless  of  where  the  spectral  peaks  are  located. 

The  assumption  of  independent  spectral  components  seems  adequate  most  of  the 
time  for  such  applications.  However,  for  very  steep  waves  or  waves  in  very  shal¬ 
low  water  the  assumption  is  often  incorrect.  Cases  have  been  documented  in  which 
all  major  spectral  components  at  frequencies  higher  than  the  main  peak  frequency 
are  closely  tied  to  the  main  peak  and  do  not  represent  independent  energy  con¬ 
centrations.  The  steepness  and  relative  depth  criteria  in  equations  (S)  and  (6) 
can  be  used  to  indicate  cases  where  the  assumption  of  independent  spectral  com¬ 
ponents  may  be  poor. 

VI .  SUMMARY 

This  report  has  presented  guidelines  for  interpreting  nondirectional  wave 
energy  spectTa,  including  a  simple  method  for  identifying  major  spectral  peaks 
and  for  estimating  significant  wave  height,  period,  and  energy  for  each  major 
peak.  Each  major  spectral  peak  is  generally  assumed  to  represent  an  independent 
wave  train.  Coarse  guidelines  are  presented  for  identifying  cases  where  major 
peaks  do  not  represent  independent  trains.  Spectral  analysis  and  display  tech¬ 
niques,  and  the  natural  variation  of  spectra  in  space  and  time,  are  discussed 
to  show  that  the  above  method  for  interpreting  spectra  provides  a  relatively 
consistent  description  of  general  spectral  characteristics.  The  method  allows 
a  more  detailed  and  accurate  description  of  ocean  surface  waves  than  that  given 
by  a  single  significant  height  and  period,  yet  it  eliminates  much  of  the  formi¬ 
dable  detail  of  a  full  spectrum. 
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